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Despite qualitative difTerences in their underlying physics, both hard and soft glassy materials 
exhibit almost identical linear rheological behaviors. We show that these nearly universal properties 
emerge naturally in a shear-transformation-zone (STZ) theory of amorphous plasticity, extended to 
include a broad distribution of internal thermal-activation barriers. The principal features of this 
barrier distribution are predicted by nonequilibrium, effective-temperature thermodynamics. Our 
theoretical loss modulus G"{uj) has a peak at the a relaxation rate, and a power law decay of the 
form cj""" for higher frequencies, in quantitative agreement with experimental data. 



Qualitatively different kinds of amorphous materials 
- e.g. structural, metallic and colloidal glasses - ex- 
hibit remarkable similarities in their linear rheological 
properties '^I despite their enormous range of internal 
dynamics and intrinsic time scales. In particular, their 
frequency dependent loss moduli G"{uj) all have peaks 
that rise near a viscous relaxation rate and drop slowly 
over many decades of higher frequencies. We show here 
that this near universality emerges naturally in the shear- 
transformation-zone (STZ) theory of amorphous plastic- 
ity , generalized to include a distribution of internal 
barriers heights of the kind found experimentally in 1980 
by Argon and Kuo f^. We also show that the principal 
features of this barrier-height distribution are predicted 
by nonequilibrium, effective-temperature thermodynam- 
ics 0. Our analysis differs from both soft glassy rheol- 
ogy (SGR) i and from mode-coupling theory [10|. In 
particular, we start with a kinematic formulation that is 
fundamentally different from that used in SGR, and we 
insist that the elementary shear tranformations be driven 
by ordinary thermal fluctuations. 

The STZ theory assumes that the degrees of freedom 
of a glassy material can be separated into two weakly 
coupled subsystems - the slow configurational degrees 
of freedom, i.e. the inherent structures, and the fast 
kinetic-vibrational degrees of freedom. While the lat- 
ter are generally in equilibrium with the thermal reser- 
voir, the former are characterized by an effective temper- 
ature that may depart from the reservoir temperature 
Q. A fundamental premise of STZ theory is that irre- 
versible shear deformations occur only at rare, localized, 
two-state, flow defects in the configurational subsystem. 
In flowing states, these STZ's appear and disappear as 
the driven system makes transitions between its inherent 
structures. In jammed states, the STZ's are configura- 
tionally frozen, but they still are able to make transi- 
tions between their internal orientations in response to 
ordinary thermal fluctuations or applied stresses. 

For spatially homogeneous systems, we may start by 
assuming for simplicity that the STZ's are oriented only 
in the ± directions relative to the shear stress s. We 
then consider STZ's characterized by an internal, thermal 



activation barrier A. Let the number of ± STZ's with 
given A be N±{A), and let the total number of (coarse- 
grained) molecular sites be N. The master equation for 
N±{A) is i 



TO 7V±(A) - R{±sc, A) N^iA) - i?(T.sc, A) iV±(A) 



p{9) [N,,{A)/2-N±{A)] 



(1) 



Here, tq is a fundamental time scale, for example, a vi- 
bration period for molecular glasses or a Brownian dif- 
fusion time for colloidal suspensions. R{±sc,A))/to is 
the rate per STZ for thermally activated transitions be- 
tween ± orientations; and sc is the partial stress acting 
on the configurational subsystem. In contrast to SGR, we 
assume that sc is a coarse-grained stress, determined by 
external forcing, and that the molecular-level stresses are 
accounted for implicitly by the two-state dynamics. The 
total stress, s = sc+'r]K * 7, includes (when appropriate) 
a partial stress acting on the kinetic-vibrational subsys- 
tem. For colloidal suspensions, this additional stress is a 
kinetic viscosity; 7 denotes the total shear rate and rjx* 
denotes a convolution over time. 

The two terms in square brackets on the right-hand- 
side of Eq. ([ij are the rates at which STZ's are cre- 
ated and annihilated by spontaneous thermal fluctua- 
tions. Mechanically generated noise, as used in [1, 0], 
is second order in the applied stress and therefore can be 
neglected in this linear theory. We are making a detailed- 
balance approximation in which Neq{A) is the value ap- 
proached by 2iV±(A) in steady-state. The sum over all 
possible A's is determined by a Boltzmann-like factor 
exp(-ez/x) according to J^a ^eq{A) = N exp{-ez /x)- 
X is the effective disorder temperature in energy units, 
and ez is a typical STZ formation energy. p(6) is the 
thermal noise strength, and 9 = kBT is the bath temper- 
ature in units of energy. 

The contribution to the rate of deformation coming 
from STZ transitions with barrier-height A is 
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eoA(A)C(sc,A) r(sc,A)-m(A) , 
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where V is the volume of the system, and vq is a molec- 
ular volume that sets the size of the plastic strain in- 
crement induced by an STZ transition. We expect cq = 
Nvq/V to be a number of the order of unity. As usual 
[dl, we have defined 



A(A) 



iV+(A) + A^_(A) 
N 



m(A) = 



iV+(A) - TV- (A) 
iV+(A) +7V_(A)' 



C(sc, A) = - [Rise, A) + Ri-sc, A) 
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(3) 



R{sc,A)+R{^sc,A)' 

According to Eq. ([ij, the equation of motion for m{A) 
is 

Tom(A) = 2C(sc, A) [r(sc, A)-m(A)] ~p{9)m{A). (4) 

In deriving Eq. Q, we have assumed that A(A)=:0. 
This means that we are focusing on time scales in which 
true structural aging is negligible, and in which the linear 
response of the system is time-translationally invariant. 
However, the onset of structural aging does play an im- 
portant role in what follows. The rate at which STZ's are 
spontaneously created and annihilated by thermal fluctu- 
ations is proportional to (p/tq) exp(— e^/x), which has 
the form of a conventional activation rate with x play- 
ing the role of the temperature and p{6)/tq being the 
attempt frequency. The kinetic prefactor p{9) is a super- 
Arrhenius function of the temperature, associated with 
the fact that the configurational rearrangements needed 
to form or annihilate STZ-likc defects involve many-body 
fluctuations that become increasingly complex and un- 
likely as the temperature decreases. Accordingly, p{9) is 
approximately equal to unity at temperatures well above 
the glass temperature, but becomes very small at lower 
temperatures, and vanishes below a glass transition tem- 
perature if such a transition exists 

We now introduce p{A), a normalized distribution over 
barrier heights such that A(A) = exp(— e2/x)p(A)dA, 
and use Eq. ^ to write the total rate of irreversible 
deformation in the form 



eoe 
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dAp{A)C{sc,A) [T{sc,A)-m{A)]. 

(5) 

To complete the derivation, we must specify the tran- 
sition rate R{sc,A), which determines C(sc,A) and 
T{sc,A) according to Eqs.([3]). In a linear theory 
C{sc, A) must be independent of sc, and T{sc, A) must 
be linear in it. The simplest assumptions for our purposes 
are 

«o sc 



Cisc,A)^Po{9,A)e-^^\ T{sc,A) 



ao I 
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Here, poid. A) is a many-body activation prefactor, which 
is analogous to p{0), but which does not necessarily van- 
ish at a glass transition, oq is a dimensionless num- 
ber of the order of unity. We stress that C{sc,A) in 



Eq. ^ describes only ordinary thermally activated pro- 
cesses. We then write i^(A) = 2po{0,A) exp(— A/6'), 
p(i^) = — p{A) dA/dv, and, in most cases, work with 
the kinetic quantity v instead of the energy A as the 
independent variable. 

To compute the linear oscillatory response, we assume 
that the total shear rate 7 is simply the sum of elastic 
and plastic parts: 



'j ^ sc/p + DP\sc), sc = s-r]K*i, 



(7) 



where p is the shear modulus. We then denote Fourier 
transforms as functions of frequency w by 7 etc., let tjk * 
7 — > i a; f}ii{ui) 7, and use the preceding equations to solve 
for G{ijj) = 5/7. The result is: 



G{lj) = IUJTq p 



iujTQ + K J{uj) 



A = e 

2006* 



°z/x 

(8) 



where 



p PTq 



J{uj) = / dvp{v) V 



ibjTo + p 



(9) 



JUJTo + p + 

The storage and loss moduli are, respectively, G' — Re[G] 
and G"' = Im[G]. 

In interpreting these results, we note first that G{ll)) in 
Eq. ([5]) cannot, in any physically meaningful way, be ex- 
pressed as an average over Maxwell modes as in SGR 
This feature is a result of our kinematic assumption in 
Eq. ([7|), plus our assumption that the plastic strain rate 
appearing there is a sum over independent contributions 
from the two-state STZ's with different A's. Next, we use 
Eq. to compute the Newtonian viscosity associated 
with configurational deformations: 
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A J(0) 



(10) 



and deduce immediately, without yet knowing p(z^), that 
the low-frequency structure of G{uj) occurs approxi- 
mately in the neighborhood of uja ^ AJ{0)/tq = p/rj^. 
Thus, Wq is the viscous relaxation rate. However, the 
structure of the a peak depends sensitively on p(z^) , and 
the approximation made above cannot replace a full eval- 
uation of G{uj). 

The distribution p{A) is a near-equilibrium property 
of the configurational subsystem; therefore, it must be 
determined by the effective temperature x- Because A 
is measured downward from some reference energy, we 
postulate, at least for some range of values of A, that 
p{A) has the form 



p(A) oc e+'^^/x. 



(11) 



where 6 is a dimensionless number. In the limit of small A 
(large z^), where po — l, Eq. PT|) becomes p(z^) ~Az^^^~'', 
where C,=hO/x and A is a normalization factor. 
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Integrability of p( A) requires that this distribution be 
cut off at some A > A*, < z/*. To estimate i/*, we 
argue that STZ transition rates cannot be slower than 
the rates at which the STZ's themselves are appearing 
and disappearing. Therefore, we propose the important 
relation: 

ly* ^2po{e,A*)e-^'/'' ^ 2 p{6) e- . (12) 

Note that we do not need to know po in order to deter- 
mine V* from the last expression on the right-hand side 

ofEq. (ini). 

Since the cutoff at A* cannot be infinitely sharp, 
we assume that the distribution drops off exponentially, 
p(A) cxexp[- (A - A*)/Ai] for A>A*, where Ai is an 
energy that may be of the order of 9. We combine these 
limiting behaviors to write 

with Ci — 9/Ai ^ 1. This is an oversimplified, three- 
parameter representation of piv) , with two exponents C, 
and Ci determining the large-z^ and small-;/ limits re- 
spectively, and a single crossover value of v* . It is re- 
markably similar to the experimentally deduced distri- 
butions shown in Q • The parameter C, controls the high- 
frequency behavior of G{u)). A simple scaling analysis 
for wtq ^ p and v* predicts that G"{uj) ^ {ujtq)^'^ 
above the a peak. 

We now compare our theoretical results to experimen- 
tal data on the rheology of both hard and soft glasses. 
We first focus on the oscillatory response of structural 
and metallic glasses, for which tq is of the order of pi- 
coseconds, wtq^I, and r]K is negligible. The interesting 
behavior occurs at temperatures near or slightly above 
the glass temperature. 

Our principal sources of information about the oscilla- 
tory responses of structural and metallic glasses are the 
papers by Gauthier et al., in particular [Ij. These au- 
thors show that the functions G{uj), for a wide variety 
of noncrystalline materials at their glass temperatures, 
have very similar behaviors. Specifically, the loss modu- 
lus G"{oj) has a broad peak at and drops off at high 
frequencies like uj~'^ as predicted above. For metallic 
glasses, Gauthier et al. find C — 0.4. 

In Fig. [H we show G'{oj)/p and G"{uj)/^ aspredicted 
by Eq. ([8]), along with data from Fig. 2 of P| for the 
metallic glass Vitreloy 4 at its glass temperature Tg. 
In estimating the theoretical parameters, we have used 
Tg ~ 600K, To ~ 2 X lO-i^sec and ~ 50 GPa. To make 
approximations for the other parameters in Eq. ([8]), we 
note that, if the volume vq is of the order of a few cu- 
bic nanometers, then the ratio vop/6g (the prefactor of 
A in Eq. (jH)) is approximately 10^. Then, to estimate 
A(x~0g) '^exp(— 62/6'^), we assume that 0g is the same 
as the steady-state value of the effective temperature. 
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FIG. 1. (Color online) Experimental data for Vitreloy 4 and 
theoretical comparisons for the storage modulus G' (uj) (red) 
and the loss modulus G"{uj) (blue). The data points were 
extracted from Fig. 2 in [Ij, where very similar curves for 
oxide and polymeric glasses can be found. 



usually denoted by xoj for systems driven persistently at 
shear rates much slower than Tq ^ . The ratio xo I ez rnay 
be a universal quantity; it usually turns out to be in the 
range 0.1 — 0.2. We therefore estimate ^g/e^^O.lS, im- 
plying that A(6'g) - 10^3, and thus that A~ 10. Then 
Eq. (Hg) tells us that v* c^XQ''^ p{e g) . 

The theoretical curves in Fig. [T]have been computed 
using Ci = 1, p{9g)/TQ = 1.25 x lO^^sec^^ and A = 25. In 
effect, we have set eo/ao ^ 2.5 in Eq. ([5]), which is well 
within our theoretical uncertainty. It is important to note 
that these parameters imply that v* ~ 10~^^, which is 
extremely small compared to its upper limit at i> = 2. So 
far as we can tell from numerical exploration, this small 
value of v* is sharply determined by the experimental 
data. The major discrepancy between the theoretical 
curve and the experimental data is that the measured 
G' {lj) is approximately linear in lj at low frequencies, 
instead of being proportional to uP as predicted by our 
theory and by Maxwell models. 

We turn now to soft glasses, in particular, to ther- 
mosensitive colloidal suspensions, whose rheology differs 
from that of structural and metallic glasses in at least 
two important respects. First, in colloidal systems, the 
approach to jamming near a glass transition is controlled 
more sensitively by the volume fraction than by the tem- 
perature. Second, the microscopically short molecular 
vibration period in structural glasses is replaced in col- 
loids by the very much longer time scale for Brownian 
motion of the particles. As a result, the high-frequency 
cutoff at 11 = 2 is probed in rheological experiments, and 
the kinetic viscosity fjK{^jj) is relevant at accessibly high 
values of uj. 

In order to evaluate G(aj) for colloidal suspensions, we 
need an expression for the frequency-dependent kinetic 
viscosity 7)x('^)- Here, we follow fl^ and write f]K{ijj) = 
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FIG. 2. Experimental data and theoretical comparisons for 
the storage modulus G' [to) (red) and the loss modulus G"(lo) 
(blue) for two different suspensions of thermosensitive parti- 
cles, as reported in [l^. The values of the parameters for the 
top and bottom panels, respectively, are: p = 0.04, 3 x 10~*; 
!/*=0.001, 10~V; C = l-0, 0.5; A = 200, 40; At=12, 35 Pa; and 
TK = 0.004, 0.002 sec. In all cases, Ci = l and c = 0.1. 



decreases slightly as the system becomes stiffer. A is very 
large for the liquidlike example in the top panel, imply- 
ing that the STZ density is large in this system. On the 
other hand, the value of A, and the relation between v* 
and p for the glassy system in the bottom panel are com- 
parable to the estimates for the bulk metallic glass in the 
preceding discussion - despite the fact that the underly- 
ing time scales for these systems differ by nine orders of 
magnitude. 

The linear oscillatory measurements discussed here are 
sensitive probes of the internal structure of glassy ma- 
terials, revealing the broad range of activation mecha- 
nisms that occur within them and the relation between 
these mechanisms and the effective-temperature thermo- 
dynamics of glassy disorder. One of the deepest questions 
that we have not addressed here, however, is the relation 
between internal STZ dynamics and true structural ag- 
ing, in which the slow configurational degrees of freedom 
relax toward thermodynamic equilibrium. We have be- 
gun to address the latter issues and will report on them 
in a forthcoming publication [14]. 

We thank M. Siebenburger for sending us the data 
shown in Fig. [2j JSL was supported in part by the 
Division of Materials Science and Engineering, Office of 
Basic Energy Sciences, Department of Energy, DE-AC05- 
OOOR-22725, through a subcontract from Oak Ridge Na- 
tional Laboratory. EB was supported in part by the 
Harold Perlman Family Foundation and in part by a 
grant from the Robert Rees Applied Research Fund. 



/iTx(c-|-ia;rx)~^^^, where tk is a viscous time scale, and 
c is a dimensionless constant that we have added in order 
to regularize this formula at small lo. 

In Fig. [2] we show two examples of how the STZ the- 
ory developed here is capable of reproducing the exper- 
imental results of Siebenburger et al. |13|. These au- 
thors explored a range of effective volume fractions 
and a wide range of frequencies w by using suspensions 
of thermosensitive particles (polystyrene cores with at- 
tached networks of thermosensitive isopropylacrylamide 
molecules). The Brownian time scale in these experi- 
ments is To ~ 0.003 sec. The volume fraction for the top 
panel is 0e// = 0.518, while that for the bottom panel 
is 0.600. The theoretical parameters, deduced by fitting 
the data, are listed in the figure caption. 

The trends are interesting. The example in the top 
panel of Fig. [2] is a system whose relatively small volume 
fraction puts it well away from the glass transition. It is 
effectively a liquid; /9 = 0.04 means that there is relatively 
little super- Arrhenius suppression of the structural relax- 
ation rate. In contrast, the system in the bottom panel 
is much more glassy, and p decreases significantly. The 
shear modulus p, increases and the viscous time scale tk 
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